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1. Introduction 

The study of event-by-event (e-by-e) fluctuations in high-energy nucleus-nucleus (A+A) col¬ 
lisions opens new possibilities to investigate properties of strongly interacting matter (see, e.g., 
Refs, [j^] and 0 and references therein). Specific fluctuations can signal the onset of deconfine¬ 
ment when the collision energy becomes sufficiently high to create the quark-gluon plasma (QGP) 
at the initial stage of A+A collision l|3|, §]. By measuring the fluctuations, one may also observe 
effects caused by the dynamical instabilities when the expanding system goes through the order 
transition line between the QGP and the hadron resonance gas Furthermore, the critical point 
(CP) of strongly interacting matter may be signaled by characteristic fluctuation pattern l|^ 0, |^. 
Therefore, e-by-e fluctuations are an important tool for the study of properties of the onset of de¬ 
confinement and the search for the CP of strongly interacting matter. However, mostly due to the 
incomplete acceptance of detectors, difficulties to control e-by-e the number of interacting nucle¬ 
ons, and also not well adapted data analysis tools, the results on e-by-e fluctuations are not yet 
mature. Even the simplest tests of statistical and dynamical models at the level of fluctuations are 
still missing. 

In this presentation the theoretical progress in several areas related to the study of e-by-e fluc¬ 
tuations is reported. A role of the global conservation laws is discussed in Sec. The strongly 
intensive measures of e-by-e fluctuations are introduced in Sec. They give a possibility to study 
e-by-e fluctuations in a physical system when its average size and size fluctuations can not be con¬ 
trolled experimentally. In Sec. ^ a novel procedure, the identity method, is described for analyzing 
fluctuations of identified hadrons under typical experimental conditions of incomplete particle iden¬ 
tification. Finally, in Sec. ^ using the van der Waals equation of state adopted to the grand canonical 
ensemble formulation we discuss particle number fluctuations in a vicinity of the CP. Most part of 
our discussion concerns the e-by-e fluctuations of hadron multiplicities. However, many of our 
physical conclusions can be applied to more general cases. 


2. Global Conservations Laws 


In this section we illustrate the role of global conservation laws in calculating of e-by-e fluc¬ 
tuations within statistical mechanics. Successful applications of the statistical model to description 
of mean hadron multiplicities in high energy collisions (see, e.g.. Refs. [^, 10, |^] and references 


therein) has stimulated investigations of properties of the statistical ensembles. Whenever possi¬ 
ble, one prefers to use the grand canonical ensemble (GCE) due to its mathematical convenience. 
The canonical ensemble (CE) should be applied [|T2, 13 ] when the number of carriers of conserved 
charges is small (of the order of 1), such as strange hadrons [ 14 ], antibaryons l|T5|], or charmed 
hadrons [p^]. The micro-canonical ensemble (MCE) has been used I ITtI , [TsI , to describe small 
systems with fixed energy, e.g., mean hadron multiplicities in proton-antiproton annihilation at rest. 
In all these cases, calculations performed in different statistical ensembles yield different results. 
This happens because the systems are ‘small’ and they are ‘far away’ from the thermodynamic 
limit (TL). The multiplicities of hadrons produced in relativistic heavy ion collisions are typically 
much larger than 1. Thus, their mean values obtained within GCE, CE, and MCE approach each 
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other. One refers here to the thermodynamical equivalence of statistical ensembles in the TL and 
uses the GCE, as a most convenient one, for calculating the hadron yields. 

A statistical system is characterized by the extensive quantities: volume V, energy E, and 
conserved charge(s)^ Q. The MCE is defined by the postulate that all micro-states with given 
V, E, and Q have equal probabilities of being realized. This is the basic postulate of statistical 
mechanics. The MCE partition function just calculates the number of microscopic states with given 
fixed (y,E,2) values. In the CE the energy exchange between the considered system and ‘infinite 
thermal bath’ is assumed. Consequently, a new parameter, temperature T, is introduced. To define 
the GCE, one makes a similar construction for the conserved charge Q: an ‘infinite chemical bath’ 
and the chemical potential /i are introduced. The CE introduces the energy fluctuations. In the 
GCE, there are additionally the charge fluctuations. Therefore, the global conservation laws of E 
and Q are treated in different ways: in the MCE both E and Q are fixed in each microscopic state, 
whereas only average value of E is fixed in the CE, and average values of E and Q in the GCE. 

The MCE, CE, and GCE are the most familiar statistical ensembles. In several textbooks (see, 
e.g.. Ref. ^), the pressure (or isobaric) ensemble has been also discussed. An ‘infinite bath 
of the fixed external pressure’ p is then introduced. This leads to volume fluctuations around the 
average value (see Ref. l^]). In general, there are 3 pairs of variables - (V^p)^ - 


and, thus, the 8 statistical ensembles^ can be constructed. 

Measurements of hadron multiplicity distributions in A+A collisions, open a new field of 
applications of the statistical models. The particle multiplicity fluctuations are usually quantified 
by a ratio of the variance to mean value, the scaled variance. 


0)[N] = 


{N^) - {Ny 

W) 


( 2 . 1 ) 


and are a subject of current experimental activities. In statistical models there is a qualitative dif¬ 
ference in properties of a mean multiplicity and a scaled variance of multiplicity distribution. It 
was recently found l^, M, 


, 33, 34] that even in the TL correspond¬ 


ing results for the scaled variance are different in different ensembles. Hence the equivalence of 
ensembles holds for mean values in the TL, but does not extend to fluctuations. Several examples 
below illustrate this statement. 


2.1 Canonical Ensemble 

Let us consider a system which consists of one sort of positively -i-l and negatively -1 charged 
particles (e.g., 7t^ and 7t~ mesons) with total charge equal to zero Q = 0 . Eor the relativistic 
Boltzmann ideal gas in the volume V at temperature T the GCE partition function reads: 

Zg,.(V,r) = t t ■ (2.2) 

iV+=0iV_=0 


^ These conserved charges are usually the net baryon number, strangeness, and electric charge. In non-relativistic 
statistical mechanics the number of particles plays the role of a conserved ‘charge’. 

^For several conserved charges {Qi} the number of possible ensembles is larger, as each charge can be treated either 
canonically or grand canonically. 
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Figure 1: Left: The ratio of (A^±)ce ( |2-6| ) to (A^±)gce as a function of z. Right: The scaled variances of 
N± calculated within the GCE, = 1 ( |2.7| ), and CE, (2.8). 


In Eq. (|2.2|) z is a single particle partition function 




/ k^dk exp 

Jo 


{k^ + m^y^'^ 


V 


T m 




(2.3) 


where m is a particle mass, and K 2 is the modified Hankel function. Parameters A+ and A_ 
are auxiliary parameters introduced in order to calculate mean numbers and fluctuations of the 
positively and negatively charged particles. They are set to one in the final formulas. The chemical 
potential equals zero to satisfy the condition (2)gee = 0. 

The CE partition function is obtained by an explicit introduction of the charge conservation 
constrain, Q = N^—N- = 0, for each microscopic state of the system and it reads: 

Zee(i/,r) = £ £ = (2.4) 

N+=0N-=0 ■'''+• 

1 


where the integral representations of the 5-Kronecker symbol and the modified Bessel function Iq 
were used. The average number of and N- can be calculated as 


(^±)gce 

(A^±)ce 




/o(2z) ■ 


(2.5) 

( 2 . 6 ) 


The exact charge conservation leads to the CE suppression, /i(2z)//o(2z) < 1, of the charged par¬ 
ticle multiplicities relative to the results for the GCE (^). The ratio of {N±) calculated in the CE 
and GCE is plotted as a function of z in Eig. left. 
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The corresponding scaled variances are [Q: 




gee 


CO" 


(jVi)gee - (A^±)^ce ^ ^ 

(A^±)gce 

{nDcc - {N±)l _ ^ \h{2z) h{2zy 

{N±)c, lIo{2z) h{2z)_ ■ 


(2.7) 

( 2 . 8 ) 


In the large volume limit (V ^ oo corresponds also to z ^ oo) one can use an asymptotic 
expansion of the modified Bessel function, 


and obtains: 


lim4(2z) 

).CX3 


exp(2z) 

\/4^ 




(2.9) 


(A^±)ce = {N±)g.c.e=Z, (2.10) 

+ = ^ CO^ . (2.11) 

The dependence of the scaled variance calculated within the CE and GCE on z is shown in Eig. D 
right. The scaled variance shows a very different behavior than the mean multiplicity: in the large 
z limit the mean multiplicity ratio approaches one and the scaled variance ratio 1/2. Thus, in the 
case of N± fluctuations the CE and GCE are not equivalent. 


2.2 Micro-Canonical Ensemble 


Our second example is the ideal gas of massless neutral Boltzmann particles. We consider the 
same volume and energy in the MCE and GCE and compare these to formulations. The fixed MCE 
energy E and the mean GCE energy (£’)gce are connected via equation (the degeneracy factor is 
^=1): 

E = (£)gee = Evt\ (2.12) 

The mean multiplicity (A/^)mce in the MCE is approximately equal to the GCE value: 

Wmee = Wgee ^ E y = N . (2.13) 


The approximation (A^)mce —N is valid for A/^ 1 and reflects the thermodynamic equivalence of 

the MCE and GCE. The scaled variances for the multiplicity fluctuations are, however, different in 


the GCE and MCE 


t^gee 

t^mce 


(A^')gce - {N)l, _ 

Wgee 

(A^^)mce - (A^)Le ^ 1 

(A^)mce 4 • 


(2.14) 

(2.15) 


Thus, despite of thermodynamic equivalence of the MCE and GCE the value of tUmce is four times 
smaller than the scaled variance of the GCE (Poisson) distribution, tOgee = 1. Eigure|^a) shows a 
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Figure 2: (a): The multiplicity distribution P{N) of massless neutral particles in the MCE, dashed line, and 
the GCE (Poisson) distribution, dotted line, for the same mean multiplicity, N = 50. (b): The momentum 
spectrum of massless neutral particles calculated within the MCE, dashed line, and the GCE (2.16), dotted 
line. The temperature is fixed as T = 160 MeV and the system energy is £" = ?>NT = 24 GeV for both plots. 


comparison of the MCE and GCE results for the multiplicity distributions P{N). Different values 
of the scaled variances ( |2.14D and ( |2.15D signify the different widths of P{N) distributions in the 
GCE abd MCE clearly seen in Fig. 

Figure ^(b) shows the single particle momentum spectra F{p) in the GCE and MCE {N — 50 
and r = 160 MeV). A single particle momentum spectrum in the GCE reads: 


Fgce{p) — 


1 dN 


V 


N p^dp 2 k^N 


f P\ 

1 / 

' P\ 

V t) 


. t) 




(2.16) 


The MCE spectrum is very close to the Boltzmann distribution ( |2.16| ) at all momenta up to 10 GeV/c 


At this momentum, the spectrum F{p) is already dropped in a comparison to its F(0) value by 
about a factor of 10“^^. The MCE spectrum decreases faster than the GCE one at high momenta. 
Close to the threshold momentum, ;? = £' = 24 GeV, where the MCE spectrum goes to zero, large 
deviations from (|2.16|) are observed. In order to demonstrate these deviations the MCE and GCE 


momentum spectra are shown in Fig. ^ (b) over 90 orders of magnitude. 


2.3 Hadron Resonance Gas in GCE, CE, and MCE 

We present now the results of the hadron resonance gas (HRG) for the e-by-e fluctuations 
of negatively charged and positively charged hadrons (see details in Ref. [@). The correspond¬ 
ing scaled variances co~ and co^ are calculated in the GCE, CE, and MCE along the chemical 
freeze-out line in central Pb-i-Pb (Au-i-Au) collisions for the whole energy range from SIS to LHC. 
The model parameters are the volume V, temperature T, baryonic chemical potential and the 
strangeness saturation parameter js- They are chosen by fitting the mean hadron multiplicities. 
Once a suitable set of the chemical freeze-out parameters is determined for each collision energy, 
the scaled variances (0~ and co~ can be calculated in different statistical ensembles. The results 
for the GCE, CE, and MCE are presented in Fig. || as functions of the center-of-mass energy 
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Figure 3: The scaled variances co~ (left) and (0+ (right), both primordial and final, along the chemical 
freeze-out line for central Pb-i-Pb (Au-i-Au) collisions. Different lines present the GCE, CE, and MCE results. 
Symbols at the curves for final particles correspond to the specific collision energies. The arrows show the 
effect of resonance decays. 




Figure 4: The scaled variance for negatively charged (left) and positively charged (right) hadrons along 
the chemical freeze-out line for central Pb-i-Pb collisions at the SPS energies. The points show the data 
of the NA49 Collaboration [|^. Curves show the GCE, CE, and MCE results calculated with the NA49 
acceptance. 

of the nucleon pair. A comparison between the data and predictions of statistical models should 
be performed for results which correspond to A-i-A collisions with a fixed number of nucleon par¬ 
ticipants. In Fig. ^ our statistical model results are compared with the NA49 data at collision 
energies 20, 30, 40, 80 and 158 A GeV for the 1% most central Pb+Pb collisions selected by the 
numbers of projectile participants. In the experimental study of A-i-A collisions at high energies 
only a fraction of all produced particles is registered. If detected particles are uncorrelated, the 
scaled variance for the accepted particles, C0acc[A^], can be presented in terms of the full An scaled 
variance as (see, e.g., & 


Wacc[^] ^ \ - q + q(o[N], (2.17) 

where ^ is a probability of a single particle to be accepted. Figure ^ demonstrates that among three 
different statistical ensembles the MCE is in a better agreement with the NA49 data. The reasons 
of this are not yet clear. 
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2.4 Generalized Statistical Ensembles 


A choice of statistical ensemble is crucial in calculating fluctuations. On the other hand, 
it is clear that GCE, CE, and MCE are only some typical examples. A general concept of the 
statistical ensembles was suggested in Ref. [^. The extensive quantities {V^E^Q) = A define the 
MCE. Different statistical ensembles are then constructed using externally given distributions of 
extensive quantities, Pa (A). The distribution of any observable O in is then obtained in two steps. 
Eirstly, the MCE 0-distribution, Pmce{0',A), is calculated at fixed values of the extensive quantities 
A. Secondly, this result is averaged over the external distribution Pa{A) 

Pa{0) = IdA Pa{A) P^eiO;A) . (2.18) 


Eluctuations of extensive quantities A around their average values depend not on the system’s phys¬ 
ical properties, but rather on external conditions. One can imagine a huge variety of these condi¬ 
tions, thus, the standard statistical ensembles discussed above are only some special examples. The 
ensemble defined by Eq. ( |2.18[ ), the a-ensemble, includes the standard statistical ensembles as the 
particular cases. The generalized statistical mechanics based on Eq. (|2.18D can be applied to dif¬ 


ferent tasks of hadron production in high energy collisions. Eor example, based on Eq. ( |2.18| ) and 
introducing the scaling volume fluctuations Pa{V), an attempt was made in Ref. l |^ ] to extend the 
statistical model to the hard domain of high transverse momenta and/or high hadron masses. 


3. Strongly Intensive Measures of E-by-E Fluctuations 

A significant increase of transverse momentum and multiplicity fluctuations is expected in the 
vicinity of the CR One can probe different regions of the phase diagram by varying the collision 
energy and the size of colliding nuclei. The possibility to observe signatures of the critical point 
inspired the energy and system size scan program of the NA61/SHINE Collaboration at the CERN 
SPS l|3^ and the low beam energy scan program of the STAR and PHENIX Collaborations at the 
BNL RHIC lpO|]. In these studies one measures and then compares e-by-e fluctuations in collisions 
of different nuclei at different collision energies. The average sizes of the created physical systems 
and their e-by-e fluctuations are expected to be rather different (see, e.g.. Ref. l|4T[]). This strongly 
affects the observed fluctuations, i.e., the measured quantities would not describe the local physical 
properties of the system but rather reflect the system size fluctuations. Eor instance, A+A collisions 
with different centralities may produce a system with approximately the same local properties (e.g., 
the same temperature and baryonic chemical potential) but with the volume changing significantly 
from interaction to interaction. Note that in high energy collisions the average volume of created 
matter and its variations from collision to collision usually cannot be controlled experimentally. 
Therefore, a suitable choice of statistical tools for the study of e-by-e fluctuations is really impor¬ 
tant. 

Intensive quantities are defined within the GCE of statistical mechanics. They depend on 
temperature and chemical potential(s), but they are independent of the system volume. Strongly 
intensive quantities introduced in Ref. I p^ are, in addition, independent of volume fluctuations. 
They are the appropriate measures for studies of e-by-e fluctuations in A+A collisions and can be 
defined for two extensive state quantities A and B. Here, we call A and B extensive when the first 
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moments of their distributions for the ensemble of possible states are proportional to volume. They 
are referred to as state quantities as they characterize the states of the considered system, e.g., final 
states of A+A collisions or micro-states of the GCE. 

There are two families of strongly intensive quantities which depend on the second and first 
moments of A and B and thus allow to study e-by-e (or state-by-state) fluctuations [^]: 


where 


A[A,B] 

^A,B] 


(5)0) [A] - (A)o)[B]J , 

{B)0)[A] + (A)0)[5] - 2((A5)-(A)(5))1 , 


0)[A] = 


{A? 


(A) 


0)5 = 


(B^) - {By 


(B) 


(3.1) 

(3.2) 

(3.3) 


and averaging (...) is performed over the ensemble of multi-particle states. The normalization 
factors Ca and Cz are required to be proportional to the first moments of any extensive quantities. 

In Ref. p3[ ] a specific choice of the Ca and Cz normalization factors was proposed. It makes the 
quantities A [A, 5] and E[A,5] dimensionless and leads to A [A, 5] = Z[A,5] = 1 in the independent 
particle model, as will be shown below. 

From the definition of A [A, 5] and F[A,5] it follows that A [A, 5] = Z[A,5] = 0 in the case of 
absence of fluctuations of A and B, i.e., for (o[A] = ft)[5] = {AB) — {A){B) — 0. Thus the proposed 
normalization of A [A, 5] and F[A,5] leads to a common scale on which the values of the fluctuation 
measures calculated for different state quantities A and B can be compared. 

There is an important difference between the F[A,5] and A[A,5] quantities. Namely, in order 
to calculate A[A,5] one needs to measure only the first two moments: (A), {B) and (A^), {B^). 
This can be done by independent measurements of the distributions Pa{^) and Pb{B). The quantity 
Z[A,5] includes the correlation term, (AB) — {A){B), and thus requires, in addition, simultaneous 
measurements of A and B in order to obtain the joint distribution Pab{A^B). 


3.1 A and Z in the Independent Particle Model 

The IPM assumes that: 

(1) the state quantities A and B can be expressed as 

A = CCi -h 0^2 +... + OCn 7 5 = ^1 “h /32 +... + Pn 7 (3.4) 

where aj and Pj denote single particle contributions to A and B, respectively, and N is the number 
of particles; 

( 2 ) inter-particle correlations are absent, i.e., the probability of any multi-particle state is the 
product of probability distributions P(a^ ,j 8 ^) of single-particle states, and these probability distri¬ 
butions are the same for all 7 = 1,..., A and independent of N, 


PN{aul5ua2,P2,---,0CN,M = ^(A^) x 5(o!i,j3i) x 5 ( 0 : 2 , 182 ) x ■ ■ ■ x P{aN,l5N) , (3.5) 


where ^{N) is an arbitrary multiplicity distribution of particles. 
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It can be shown p3[ ] that within the IPM the average values of the first and second moments 
of A and B are equal to: 


(A) = a (A), 

(A^) = a2 [(A^) 

- {N)] , 

(3.6) 

(B) = P{N), 

(tf ) =W{N)+ f [{W") 

- {N)] , 

(3.7) 

(AB) = ap {N) + a- p [(A^) - (A)] . 


(3.8) 


The values of (A) and {B) are proportional to the average number of particles {N) and, thus, to the 
average size of the system. These quantities are extensive. The quantities a, j8 and j8^, aj8 
are the first and second moments of the single-particle distribution P(a,j8). Within the IPM they 
are independent of {N) and play the role of intensive quantities. 

Using Eq. ( ^ ) the scaled variance co[A] which describes the state-by-state fluctuations of A 
can be expressed as: 


{A2) _ (A)2 

^ (A) 


a2 - 


_ / A ^ 2 \ _ /^\2 _ 

=- ^ - - /AJ\ ^ ^ ’ (3.9) 

a (N) 

where a)[a] is the scaled variance of the single-particle quantity a, and 0)[N] is the scaled variance 
of N. A similar expression follows from Eq. ( p7^ for the scaled variance (o[B]. The scaled variances 
a)[A] and a)[6] depend on the fluctuations of the particle number via a)[A]. Therefore, a)[A] and 
(o[B] are not strongly intensive quantities. 

From Eqs. one obtains expressions for A[A,fi] and E[A,fi], namely: 


A[A,B] = /3 (o[a] - a a)[j3] 

Ca L 


L[A,B] = 

Thus, the requirement that 


M 

Ce 


/3 (o[a] -I- a (o[j3] - 2 ( a/3 - a • /3 


(3.10) 

(3.11) 


within the IPM leads to: 

Ca 

Ce 


A[A,B] = L[A,B] = 1 , 


(A) 

(A) 


P (o[a] - a 0)[P] , 
j8 ft)[o:] -f a (o[P] - 2 


ap 


a-p 


(3.12) 

(3.13) 

(3.14) 


In the IPM the A and B quantities are expressed in terms of sums of the single particle vari¬ 
ables, a and p. Thus in order to calculate the normalization Ca and Ce factors one has to mea¬ 
sure the single particle quantities a and /3. However, this may not always be possible within a 
given experimental set-up. For example, A and B may be energies of particles measured by two 
calorimeters. Then one can study fluctuations in terms of A [A, 6] and E[A,6] but can not calculate 
the normalization factors which are proposed above. 

We consider now two examples of specific pairs of extensive variables A and B. In the first 
example, we use the transverse momentum Pj — +. ■■pf\ where p^'* is the absolute value of 

the i* particle transverse momentum, and the number of particles A. The requirement that 


A[Pt,N] = E[Pr,A] = 1 


(3.15) 
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for the IPM leads then to the normalization factors 


[43] 


Ca = Ce = (o[pt\ ■ {N) , co[pt] = 


(3.16) 


where (o\pt\ describes the single particle /^^-fluctuations. 

As the second example let us consider the multiplicity fluctuations. Here A and B will denote 
the multiplicities of hadrons of types A and B, respectively (e.g., kaons and pions). One obtains 

ii 


Ca^{B)-{A), Cz^{A) + {B) 


(3.17) 


The normalization factors ( 3.16D and ( |3.17D are suggested to be used for the calculation for A 


and Z both in theoretical models and for the analysis of experimental data (see Ref. l |43n for further 
details of the normalization procedure). 

The <I> measure, introduced some time ago Q], belongs to the Z family within the current 
classification scheme. The fluctuation measure was introduced for the study of transverse mo¬ 
mentum fluctuations. In the general case, when A = X represents any motional variable and B — N 
is the particle multiplicity, one gets: 


= 


xco\x\ 


1 V2 r 


vnx,N] - 1 


(3.18) 


For the multiplicity fluctuations of hadrons belonging to non-overlapping types A and B the con¬ 
nection between the ^[A,5] and Z[A,5] measures reads: 


^[A,B] = 


(A) + {B) 


vnA,B] - 1 


(3.19) 


The IPM plays an important role as the reference model. The deviations of real data from 
the IPM results Eq. ( |3.15| ) can be used to learn about the physical properties of the system. This 
resembles the situation in studies of particle multiplicity fluctuations. In this case, one uses the 
Poisson distribution P{N) = exp (—X) /N\ with ft)[A] = 1 as the reference model. The other 
reference value co[N] =0 corresponds to A = const, i.e., the absence of A-fluctuations. Values of 
ft)[A] > 1 (or ft) [A] » 1) correspond to “large” (or “very large”) fluctuations of A, and ft)[A] < 1 
(or ft) [A] <C 1) to “small” (or “very small”) fluctuations. 

The fluctuation measures A and Z do not depend on the average size of the system and its 
fluctuations in several different model approaches, namely, statistical mechanics within the GCE, 


model of Multiple Independent Sources, Mixed Event Model (see Ref. l |43p for details). Thus, 
one may expect that A and Z preserve their strongly intensive properties in many real experiments 
too. An extension of the strongly intensive fluctuation measures for higher order cumulants were 
suggested in recent paper l|45[]. 

3.2 A and Z Evaluated in Specific Models 

The UrQMD l^] calculations of A and Z measures were performed in Refs. The 

Monte Carlo simulations and analytical model results for A[P 7 ’,iV] and Z[TV,iV] were presented 
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in Ref. p9|]. These measures were also studied in Ref. [ pq ] for the ideal Bose and Fermi gases 
within the GCE. The GCE for the Boltzmann approximation satisfies the conditions of the IPM, 


i.e., Eq. ( |3.12| ) is valid. The following general relations were found 




(3.20) 

(3.21) 


i.e., the Bose statistics makes A[P 7 ’,iV] smaller and r[P 7 ’,iV] larger than unity, whereas the Fermi 
statistics works in the opposite way. The Bose statistics of pions appears to be the main source 
of quantum statistics effects in a hadron gas with a temperature typical for the hadron system 
created in A+A collisions. It gives about 20% decrease of A[P 7 , A] and 10% increase of A], 


at r = 150 MeV with respect to the IPM results ( |3.12| ). The Fermi statistics of protons modifies 
insignificantly A[Pr,A] and Z[Pr,A] for typical values of T and jis. Note that UrQMD takes 
into account several sources of fluctuations and correlations, e.g., the exact conservation laws and 
resonance decays. On the other hand, it does not include the effects of Bose and Fermi statistics. 
First experimental results on A[Pr, A] and A] in p+p and Pb+Pb collisions have been reported 
in Refs. 1^. 

Now we consider the A and Z measures for two particle multiplicities Ai and A 2 . Resonance 
decays, when particle species 1 and 2 appear simultaneously among the decay products, lead to the 
(positive) correlations between Ai and A 2 numbers. Let Ai = and A 2 = 7r~ are the multiplicities 
of positively and negatively charged pions, respectively. A presence of two components is assumed: 
the correlated pion pairs coming from decays, 7r+ + 7r“, and the uncorrelated 7r+ and 7t~ 

from other sources. The 7r+ and numbers are then equal to: 


7r~^ — n-\- -|- R 


7171 5 


TT = ti- + R 


7171 5 


(3.22) 


where R^^ is the number of resonances decaying into 7i^7t~ pairs, while ^ 2 + and n- are the numbers 
of uncorrelated 7r+ and 7r“, respectively. 

Using approximate relations. 


(0[k'^] = (0[k ] = (o[Rnn] = 1 , 


one obtains two alternative expressions for the number of R;i: 7 r-i'esonances [|54|] : 


{Rn 


{k-) + {7Z+) 

{Run) 

{n-) + {K+) 


{K+Tl ) - {k+){tI ) _ 

—{FTmF)— = Pi’' 1 

1 - Z[;r+,7r-] 


(3.23) 

(3.24) 

(3.25) 


i.e., R;i: 7 r can be calculated using the measurable quantities of e-by-e fluctuations. As a result, these 
resonance abundances, which are difficult to be measured by other methods, can be estimated by 
measuring the fluctuations and correlations of the numbers of stable hadrons. Note that an idea to 
use the e-by-e fluctuations of particle number ratios to estimate the number of hadronic resonances 


was suggested for the first time in Ref. l|55[]. 

The analysis of and k~ fluctuations and correlations from resonance decays is done in 


Ref. 154]. It is based on the HRG model, both in the GCE and CE, and on the relativistic transport 
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model UrQMD. This analysis illustrates a role of the centrality selection, limited acceptance, and 
global charge conservation in A+A collisions. We present now the results for Eq. (|3.24| ) and ( |3.25|) 
in Pb+Pb and p+p collisions. Note that p in Eq. ( |3.24D is an intensive but not strongly intensive 
quantity. Thus, it is expected to be sensitive to the system size fluctuations. On the other hand, the 
Z in Eq. ( |3.25| ) is the strongly intensive measure and it should keep the same value in a presence of 
the system size fluctuations. 

The samples of 5% central Pb+Pb collision events at ^/sm = 6.27 and 17.3 GeV are consid¬ 
ered. These UrQMD results are compared with those for most central Pb+Pb collisions at zero 
impact parameter, Z? = 0 fm, and for p + p reactions at the same collision energies. Several mid¬ 
rapidity windows — Ay/2 < y < Ay/2 for final 7r+ and 7t~ particles are considered. 

A width of the rapidity window Ay is an important parameter. The two hadrons which are 
the products of a resonance decay have, in average, a rapidity difference of the order of unity. 
Therefore, while searching for the effects of resonance decays one should choose Ay > 1 to enlarge 
a probability for simultaneous hit into the rapidity window Ay of both correlated hadrons (e.g., 
7r+ and 7t~) from resonance decays. Thus, Ay should be large enough. However, Ay should be 
small in comparison to the whole rapidity interval AT /m) accessible for final hadron 

with mass m. Only for Ay <C AT one can expect a validity of the GCE results. Considering a 
small part of the statistical system, one does not need to impose the restrictions of the exact global 
charge conservations: the GCE which only regulates the average values of the conserved charges 
is fully acceptable. Eor large Ay, when the detected hadrons correspond to an essential part of the 
whole system, the effects of the global charge conservation become more important. In the HRG 
this should be treated within the CE, where the conserved charges are fixed for all microscopic 
states. The global charge conservation influences the particle number fluctuations and introduces 
additional correlations between numbers of different particle species. 

The UrQMD values of correlation parameter p[7r+, 7t~] in 5% central Pb+Pb collision events 
are shown by full circles in Pig. ||(a) and full boxes in §(b) as functions of the acceptance windows 
Ay. The UrQMD results for the most central Pb+Pb collision events with zero impact parameter, 
Z? = 0 fm, are shown by open symbols. The triangles show the results of the UrQMD simulations 
'm p-\- p reactions. The collision energy is taken as — 6.27 GeV in Pig. ^(a) and — 

17.3 GeV in Pig. |^b). The windows at the center of mass mid-rapidity are taken as Ay = 0.2, 
1, 2, 3, 4, and oo. A symbol oo denotes the case when all final state particles are detected (i.e., a 
full 4;r-acceptance). Note that the UrQMD model does not assume any, even local, thermal and/or 
chemical equilibration. Therefore, a connection between the UrQMD and HRG results for particle 
number fluctuations and correlations is a priori unknown. 

Eor the 5% most central Pb+Pb collisions the correlation p[7r+,7r“] increases with Ay. As 
seen from Pig. ||(b), this increase is rather strong at high collision energy: at large Ay, the value of 
p[7r+, 7t~] becomes much larger than the HGM results in both the CE and GCE. The behavior of 
and (o[k~] is rather similar to that of p[7r+, K~]. 

Selecting within the UrQMD simulations the most central Pb+Pb collision events with zero 
impact parameter Z? = 0 fm, one finds essentially smaller values of p[;r+,;r“] (open symbols in 
Pig. This means that in the 5% centrality bin of Pb+Pb collision events large fluctuations of the 
number of nucleon participants (i.e., the volume fluctuations) are present. These volume fluctua¬ 
tions produce large additional contributions to the scaled variances of pions and to the correlation 
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Figure 5: The UrQMD results for p[7r+, 7i~] in central Pb+Pb and inelastic p-\- p collisions for the mid¬ 
rapidity windows Ay, i.e., the center of mass rapidities y of final and satisfy the condition —AT < 
y < AT. Full circles and squares correspond to the 5% centrality selection and the open ones to the most 
central Pb-^Pb collision events with zero impact parameter Z? = 0 fm. Triangles correspond to the UrQMD 
simulations of inelastic p^p interactions. The collision energies are: s/s^ = 6.27 GeV in {a) and s/s^n = 
17.3 GeV in (b). The horizontal dashed and dashed-dotted lines show, respectively, the GCE and CE results 
taken at the corresponding ^/sm• 

parameter p[7r+,7r“]. They become more and more important with increasing collision energy. 
This is due to an increase of the number of pions per participating nucleon with increasing col¬ 
lision energy. However, one hopes that these volume fluctuations will be canceled out to a large 
extent when they are combined in the strongly intensive measures. Note also that the UrQMD 
results for p[;r+,;r“] in inelastic p + p collisions, shown in Fig. ^by triangles, are qualitatively 
similar to those in Pb+Pb collisions at Z? = 0 fm. 

The UrQMD results for E[;r+,7r“] in Pb+Pb collisions at ^/snn = 6.27 GeV and 17.3 GeV 
are presented in Fig. ^(a) and (b), respectively, as a function of the acceptance window Ay at mid¬ 
rapidity. In contrast to the results shown in Fig. both centrality selections in Pb+Pb collisions 
(5% centrality bin and Z? = 0 fm) lead to very similar results for E[7r+, 7t~] shown in Fig. ^ This 
means that the measure Z[7r+,7r“] has the strongly intensive properties, at least in the UrQMD 
simulations. The UrQMD results in p+p reactions are close to those in Pb+Pb ones. 

The GCE and CE results of the HRG are presented in Figs. ^ and by the horizontal dashed 
and dashed-dotted lines, respectively. The UrQMD results for Z[7r+,7r“], presented in Fig. 
demonstrate a strong dependence on the size of rapidity window Ay. At Ay = 1 these results are 
close to those of the GCE HRM. On the other hand, with increasing Ay the role of exact charge 
conservation becomes more and more important. From Fig. one observes that the UrQMD values 
of E[;r+, 7t~] at large Ay are close to the CE results. 

As seen from Fig. a similar correspondence between the UrQMD results for p[7r+, 7r~] in 
Pb+Pb collisions at Z? = 0 fm and their GCE and CE values is approximately valid. However, this 
is not the case for the 5% most central Pb+Pb events. In that centrality bin the volume fluctuations 
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Figure 6: The same as in Fig. ||but for K ]. 


give the dominant contributions to p [;r+, K ] for large Ay. 

For very small acceptance, Ay <C 1, one expects an approximate validity of the Poisson dis¬ 
tribution for any type of the detected particles. Their scaled variances are then close to unity, i.e., 
(o[k~] = = 1. Particle number correlations, due to both the resonance decays and the global 

charge conservation, become negligible, i.e., p[;r+,7r“] <C 1. These expectations are indeed sup¬ 
ported by the UrQMD results at Ay = 0.2 presented in Fig. Therefore, one expects Z[7r+, K~] 1 

at Ay ^ 0. This expectation is also valid, as seen from the UrQMD results at Ay = 0.2 presented in 
Fig.g 


4. Incomplete Particle Identification 

Fluctuations of the chemical (particle-type) composition of hadronic final states in A+A col¬ 
lisions are expected to be sensitive to the phase transition between hadronic and partonic matter. 
First experimental results on e-by-e chemical fluctuations have already been published from the 
CERN SPS and BNL RHIC, and more systematic measurements are in progress. 

Studies of chemical fluctuations in general require to determine the number of particles of 
different hadron species (e.g., pions, kaons, and protons) e-by-e. A serious experimental problem in 
such measurements is incomplete particle identification, i.e., the impossibility to identify uniquely 
the type of each detected particle. The effect of particle misidentification distorts the measured 
fluctuation quantities. For this reason the analysis of chemical fluctuations is usually performed in 
a small acceptance, where particle identification is relatively reliable. However, an important part 
of the information on e-by-e fluctuations in full phase space is then lost. 

Although it is usually impossible to identify each detected particle, one can nevertheless de¬ 
termine with high accuracy the average multiplicities (averaged over many events) for different 
hadron species. 
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4.1 The Identity Method 


The identity variable was introduced in Ref. []56[], and in Ref. [ p7| ] a new experimental tech¬ 
nique called the identity method was proposed. It solved the misidentification problem for one 
specific combination of the second moments in a system of two hadron species (‘kaons’ and ‘pi- 
ons’). In Refs. [^,|5^] this method was extended to show that all the second moments as well as the 
higher moments of the joint multiplicity distribution of particles of different types can be uniquely 

can 


reconstructed in spite of the effects of incomplete identification. Notably, the results [|^, | 
be used for an arbitrary number ^ > 2 of hadron species. It is assumed that particle identifica¬ 
tion is achieved by measuring the particle mass m. Since any measurement is of finite resolution, 
we deal with continuous distributions of observed masses denoted as pj{m) and normalized as 

a=l,...,/:>2) 


/ 


dmpj{m) — {Nj 


(4.1) 


Note that for experimental data the functions pj{m) for particles of type j are obtained from the in¬ 
clusive distribution of the m-values for all particles from all collision events. The identity variables 
wj(m) are defined as 


Wj{m) = 


Pjjm) 

p{m) 


pH = LP'H 

i—\ 


(4.2) 


Complete identification (Cl) of particles corresponds to distributions pj{m) which do not overlap. 
In this case, Wj = 0 for all particle species i ^ j and wj = 1 for the jth species. When the distri¬ 
butions pj{m) overlap, Wj{m) can take the value of any real number from [0,1]. We introduce the 
quantities 




N{n) 

L Wjinti) , 
HI 


wf 



WpWq 



with j — 1, • • • , ^ and I < p < q <k, and define their event averages as 


I New 1 

> (w) = ^Lw> 




N, 


ev n=l 


Ne 


ev n=l 


(4.4) 


where N^y is the number of events, and N{n) = Ni{n)^ - \-Nk{n) is the total multiplicity in the 

nih event. Each experimental event is characterized by a set of particle masses {mi,m 2 ,... ^m^}, 
for which one can calculate the full set of identity variables: {wy(mi), Wy (m 2 ),..., Wy(miv)}, with 
j = 1,... ,^. Thus, the quantities Wj, Wj, and WpWq are completely defined for each event, and 
their average values can be found experimentally by straightforward e-by-e averaging. The 
main idea is to find the relations between these W-quantities and the unknown moments of the 
multiplicity distribution {N'j) and {NpNq). In the case of Cl, one finds Wj = Nj, thus, Eq. 
yields 


{Wf) = {Nj) , {WpWj) = {NpNj) . 


(4.5) 
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4.2 Second Moments of Chemical Fluctuations 


The quantities (Wj) and {WqWp) can be calculated as follows \ p^ 


k k 

{Wf) ^ - {ujy]+'^{Nf)iujy + 2 £ {NiN,)ujiUj, , 

i=l i=l l<i<l<k 


(WpW,) = £(A^,) 




+ Zi^bupiUki + £ (NiNi) 

i—l l<i<l<k 




(4.6) 


(4.7) 


In Eqs. ( |4^ and .. ,Nk) is the multiplicity distribution, Pi{m) = pi{m)/{Ni) are the 

mass probability distributions of the ith species, and (s — 1,2) 


)i = 7^ Jdmwym) pi{m) , Upqi = 7^ Jdmwp{m)wq{m) pi{m) . 


(4.8) 


In the case of Cl, when the distributions pj{m) do not overlap, one finds that 

— ^ji 5 ^pqi — 0 5 (4.9) 

and Eqs. and ^7^ reduce then to Eq. The incomplete particle identification transforms 
the second moments (Nj) and (NpNq) to the quantities (IT?) and (WpWq), respectively. Each of the 
later quantities contains linear combinations of all the first and second moments, (Nt) and (Nf), as 
well as all the correlation terms {NiNi). Having introduced the notations 


(Wf) - [4 - [ujif] = bj , {WpWq) - Y^iNi) [, 

i—l i—l 

one can transform Eqs. and ( |47^ to the following form: 


) Uji + 2 ^ [NiNi) UjiUji — bj ^ j — 1,2,... ,A:, 

k 

I 

1=1 


^pqi ^pi^qi\ — ’ 


(4.10) 

(4.11) 

(4.12) 


) ^pi^qi~^ ^ i^i^l) i^pi^ql ^pl^qij — ^pq •> \ ^ p q ^ k . 

1=1 l<i<l<k ^ ^ 

The right-hand side of Eqs. ( |4.1lD and ( |4.12| ) defined by Eq. ( |4.10| ) are experimentally measurable 

quantities. The same is true for the coefficients - (with ^ = 1 and 2) entering the left-hand side of 

Eqs. ( |4.1lD and ( K.12D . Therefore, Eqs. ( K.HP and ( |4.12| ) represent a system of k-^k{k— l)/2 linear 
equations for the k second moments {N^) with j = 1,... and k{k — l)/2 correlators {NpNq) with 
\ < p < q <k. In order to solve Eqs. ( |4.1ip and ( 4.12 ) we introduce the [k + k{k— l)/2] x [k + 
^(^ — 1) /2] matrix A 


( 


A = 


n2 


n2 


n2 


"(i^-i)^ I ^\k-i)k 




,12 




(k-\)k 


62 


p-^)k 

^12 


(k-i)k 

^(k-l)k/ 


(4.13) 
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where 


= ^ji 

a'p^ = UniU. 

Am _ 


1 <ij <k\ 


Cli — ^UipUiq , 


\ < p < q <k ^ ; 


pi^qi 5 


I < p < q < k , i = I,... ,k ; 


pq — ^pl^qm ^ql^pm i P^q — k^ \ ^ I ^ fTl ^ k . 


(4.14) 

(4.15) 

(4.16) 


The solution of Eqs. ( |4.11D and ( 4.12 ) can be presented by Cramer’s formulas in terms of the 
deter m inants 


, detA; 


{NM = 


detA 


pq 


detA 


(4.17) 


where the matrices Aj and Apq are obtained by substituting in the matrix A the column aj,..., , a\ 2 , 

• • • ^^lk-i)k column ,... respectively, for the column b\^...^bk^ b^^ 

• • • Therefore, if detA ^ 0, the system of linear equations ( |4.11[ ) and ( |4.12D has a unique 

solution ( 4-.17 ) for all the second moments. In the case of Cl (^), one finds detA = 1, detA^ = bj, 
and detAp^ = bpq. The solution (|4.17| ) reduces then to Eq. ( |43| ). 

Introducing the [k + k{k — l)/2]-vectors 


^ = 


( (^ 1 ) \ 


( bx \ 

H) 

, ^ = 

bk 

{N 1 N 2 ) 


b\2 

\{Nk-iNk)) 


\b(k-i)k) 


(4.18) 


one can write Eqs. ( |4.11[ ) and ( 4-.12D in the matrix form A^ = The solution ( 4-.17 ) can be then 
rewritten as 


^ = A-i^, 


(4.19) 


where A ^ is the inverse matrix of A. Eor two particle species, k —2, this solution takes the form 


(^i) \ 

{nD 

\{NiN2) 


Then Eq. ( |4.20[ ) yields 


( U\y 

“21 

\MiiM21 


“12 

U 22 

U\2U22 


2u\\U\2 ^ 

2u2\U22 

U\lU22 + U\2U2l J 


-1 


(h\ 

bi 

\bn) 


{Nh 

{nD 


{N 1 N 2 ) 


b\U22 + ^2Mi2 - '^b\2U\2U22 

{u\\U 22 — Ul 2 U 2 \)^ 

Z 72 M 11 + ^1M21 “ 2Z712M21M11 
(m 11M22 — Ml2M2l)^ 

^12(m11M22 + M12M2i) -^iM22M21 -^2MiiMi2 
(m11M22 — M12M2i)^ 


(4.20) 


(4.21) 

(4.22) 

(4.23) 
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The above procedure eliminates the effect of misidentification and provides the values of all the sec¬ 
ond moments {Nj) and {NpNq) in a model-independent way, as they would be obtained in an exper¬ 
iment in which each particle is uniquely identified. This method was generalized to determine third 
and higher moments of the multiplicity distributions in events consisting of an arbitrary number of 
different particle species in Ref. [|5^. Measurements of the third and higher moments of e-by-e 
fluctuations are expected to be more sensitive for the search of the CP in A+A collisions @01. 
First results on fluctuations based on the identity method were presented in Refs, l^, |^, 


5. Fluctuations at the Critical Point 


In this section we consider the van der Waals (VDW) equation of state with both repulsive 
(b > 0) and attractive (a > 0) terms (see, e.g.. Refs. l |6^ , |65|]): 


p(yj,N) 


NT 

V - bN 


nT 

I — bn 


an 


(5.1) 


It gives an example of a system with the CP. The equation of state ( ^T|) was suggested in 1873, 
and for his work van der Waals obtained the Nobel Prize in physics in 1910. The particle proper 
volume parameter in Eq. ( |5T| ) equals iob — A- {Anr^/3) with r being the corresponding hard sphere 
radius of particle. The b parameter describes the repulsion between particles and can be rigorously 
obtained (in particular, a factor of 4 in the expression for b) for a gas of the hard balls at low density 
(see, e.g., Ref. @). 


5.1 Phase Diagram 


The VDW equation of state contains a order liquid-gas phase transition and CP. The ther¬ 
modynamical quantities at the CP are equal to [1^]: 


8a 


1 


21b ’ 3(7 ’ 27(72 ■ 


(5.2) 


The VDW equation ( |5T| ) can be then rewritten in the following dimensionless (reduced) form: 

^fn 


3 — n 


— , 


(5.3) 


where n — n/n^ p — p/Pc, and T — T/Tc. In the dimensionless presentation ( p^ the VDW 


equation has a universal form independent of the values of a and b, and the CP ( p.2D is transformed 

to fc = Pc = = 1 • 

The dimensionless VDW isotherms are presented in Fig. ^ {left) as functions ofv = n~^ . To 
describe the phase coexistence region below the critical temperature the VDW isotherms should 
be corrected by the well known Maxwell construction of equal areas. These corrected parts of the 
VDW isotherms are shown by the solid horizontal lines. Figure ^ {right) depicts the liquid-gas 
coexistence region on the in^T) plane. 
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Figure 7: Left. The dimensionless form of the VDW isotherms for pressure, p, versus the volume per 
particle v = n~^. The dashed-dotted lines present the metastable parts of the VDW isotherms at T < 1, 
whereas the dotted lines correspond to unstable parts. The full circle on the T = 1 isotherm corresponds to 
the CP. Right The phase diagram for the VDW equation of state on the in^T) plane. The phase coexistence 
region resulting from the Maxwell construction is depicted by grey shaded area. 


5.2 The GCE Formulation and Particle Number Fluctuations 


As the first step, the VDW equation (|5T| ) will be transformed to the GCE. The GCE for¬ 
mulation will be then used to calculate the particle number fluctuations in the VDW gas. These 
calculations were done in Ref. [1^, where details can be found. The GCE form of the VDW 
equation of state has its simplest form for the particle number density, 


n(7’,ju) 


Wid(r,At*) 

1 + bniii{T, n*) 


li* ^ II - b 


nT 

I—bn 


+ 2an , 


(5.4) 


as a function of temperature T and chemical potential ji. If a = Z? = 0, the function ( ^ ) is reduced 
to the ideal gas particle number density. Note that ji regulates {N)gcQ, i-e., N plays the role of 
conserved charge, and no antiparticles are introduced. Eor the ideal Boltzmann gas one obtains 





(5.5) 


where g is the degeneracy factor, m is the particle mass, and K 2 is the modified Hankel function. 
The function is the solution of transcendental equation ( |5^ . The GCE pressure 

for the VDW equation of state is obtained substituting n in Eq. ( ^d] ) by function n(T,iu) ( |5^ . 

The variance of the total particle number in the GCE can be calculated as 


Var[N] = (N^) 




(5.6) 


where symbol (...) denotes the GCE averaging, and n{T,^) is the particle number density in the 
GCE. The scaled variance for the particle number fluctuations is then |^]: 

, , VarlN] T (dn\ \ 1 2an 

It is clearly seen from Eq. ( |577| ) that in the VDW gas the repulsive interactions suppress the particle 
number fluctuations, whereas the attractive interactions lead to their enhancement. Note that for 


(5.7) 
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Figure 8: Left. The lines of constant values of the scaled variance co[N] are shown on the (n,r) phase 
diagram, outside of the mixed phase region. Right The lines of constant values of the scaled variance co[N] 
are shown on the (n, T) phase diagram for both stable and metastable pure phases. The boundary between 
stable and metastable phases is depicted by the thick black line, and the unstable region is depicted by the 
gray area. 


a = 0 the scaled variance ( |5^ is reduced to the result for the excluded volume model obtained 
earlier in Ref. [j^. 

The scaled variance ( |5^ expressed in terms n and T equals to 

1 -1 


m[7V] = 


1 


1 


9 [(3 - 


n 

AT 


(5.8) 


In Fig. P the lines of constant values of co[N] are shown on the (n^T) phase diagram outside of 
the mixed phase region. At any fixed value of T, the particle number fluctuations ( ^ ) approach 
to those of the ideal gas, i.e., ft)[A] = 1, at n ^ 0, and become small, ft)[A] <C 1, at n ^ 3. As it 
should be, the scaled variance is always positive for all possible values of n and T outside of 
the mixed phase region. At T ^ 0 this is insured by simultaneous limiting behavior of n ^ 0 in the 
pure gaseous phase and n ^ 3 in the pure liquid phase. 


5.3 Critical Point 


At the CP (r = n = 1) the scaled variance of the particle number fluctuations diverges in the 
GCE. To study the behavior of ft) [A] in a vicinity of the critical point we introduce the quantities 
T = r — 1 <C 1 and p = n — 1 <C 1. Expanding ( ^ ) at small T and p and keeping only the lowest 
orders of them one finds 


co[N] ^ 


^ + 4P + 


-1 


In particular, 


ft)[A] = ^ at p=0, and ft)[A] = ^p ^ at T = 0 

9 27 


(5.9) 


(5.10) 


Note that thermodynamical parameters T and n correspond to points outside the mixed phase re¬ 
gion. This is shown in Fig. ^ left. Thus, in Eqs. (|5^ and ( |5.10| ), parameter T is positive, while p 
can be both positive and negative. 
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Figure 9: The lines of constant values of co[N] are shown on the phase diagram for 0 < f <3 and for all 
possible n values, for stable states both inside and outside the mixed phase region. 



The VDW equation of state permits the existence of metastable phases of super-heated liquid 
and super-cooled gas. These states are depicted by the dash-dotted lines on the VDW isotherms in 
Fig. ^ {left). In metastable phases the system is assumed to be uniform and, therefore, one can use 
Eq. to calculate particle number fluctuations in these phases. 

The lines of constant values of co[N] are shown on the {n^T) phase diagram in Fig. | right 
for both stable and metastable pure phases, while the unstable region is depicted by the gray area. 
It is seen that the scaled variance remains finite, and diverges only at the boundary between the 
metastable and unstable regions. We recall that at this boundary dpjdn^ 0, where p is the di¬ 
mensionless CE pressure (5.3). One can easily show using Eqs. ( |0[ ) and (5.8) that ft)[V] ^ oo 
when dpjdn — 0. Note that metastable regions of the equation of state can be reached within fast 
non-equilibrium processes, whereas the unstable region is physically forbidden. Note also that the 
thermodynamical relations are not fulfilled in the unstable region, e.g., nonphysical behavior with 
co[N] <0 is found in this region. 

Finally, the scaled variance co[N] is shown in Fig. ^ for all possible values of n and T, both 
inside and outside of the mixed phase. Now only thermodynamically stable states are considered. 
Details of the calculations inside the mixed phase are presented in Ref. [|6^. In a vicinity of the 
critical point inside the mixed phase we introduce 0<^ = 1 — T <Cl and find at ^ ^ 0, 


co[N] ^ 


16 


.-1 


(5.11) 


Thus, the scaled variance co[N] diverges at the critical point reached from both outside and inside 
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the mixed phase. 

6. Summary 

1. The global conservation of energy and conserved charges influences the measured e-by-e 
fluctuations of all observables. Thus, to calculate these fluctuations within statistical mechanics one 
should properly choose the statistical ensemble. The GCE, CE, and MCE are only some particular 
examples. Real situation may correspond to very different externally given distributions of the 
volume y, energy E, and conserved charge(s) Q. 

2. To avoid the trivial contributions from the system size fluctuations the strongly intensive 
measures of e-by-e fluctuations should be used. These A [A, 5] and E[A,5] measures are constructed 
from the second moments (A^), and (AB) of two extensive quantities A and B. 

3. A special procedure, the identity method, should be used to calculate the chemical fluctua¬ 
tions in a case of the incomplete particle identification. It solves the misidentification problem for 
the second and higher moments. Thus the joint multiplicity distribution of particles of different 
types (a number of species k>2) can be uniquely reconstructed. 

4. The van der Waals equation of state is considered in the GCE, and particle number fluc¬ 
tuations are calculated. The scaled variance ft) [A] diverges at the the critical point. Admitting 
a presence of metastable states, one observes that the scaled variance ft) [A] diverges also at the 
boundary between the metastable and unstable regions, where dp jdn — Q. 
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